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Abstract 

The electron transport through the nanotube junctions which connect the 
different metallic nanotubes by a pair of a pentagonal defect and a heptagonal 
defect is investigated by Landauer's formula and the effective mass approx- 
imation. From our previous calculations based on the tight binding model, 
it has been known that the conductance is determined almost only by two 
parameters, i.e., the energy in the unit of the onset energy of more than two 
channels and the ratio of the radii of the two nanotubes. The conductance is 
calculated again by the effective mass theory in this paper and a simple an- 
alytical form of the conductance is obtained considering a special boundary 
conditions of the envelop wavefunctions. The two scaling parameters appear 
naturally in this treatment. The results by this formula coincide fairly well 
with those of the tight binding model. The physical origin of the scaling law 

is clarified by this approach. 
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Since the discovery of the carbon nanotubes [0] , they attract much interest as one dimen- 
sional conductor with nanometer size. It has been theoretically predicted that they become 
metallic or semiconducting according to the radius and the helicity of the honeycomb lattice 
forming the tubes. [Q] Electronic devices with nanometer size might be available by using 
the metallic nanotubes as leads. Therefore we concentrate our discussion to the metallic 
nanotubes in this paper. Recently, it becomes possible to measure conductance of the indi- 
vidual nanotubes. |§ It is expected that such technological progress leads the realization of 
the nanometer devices. 

The nanotubes can be combined naturally by introducing a pair of a pentagonal defect 
and a heptagonal defect as observed by Iijima. Such defects are called disclinations 

and are necessitated to form various structures composed of curved surface of graphitic 
layer, for example, fullerenes, || minimal- surface structures torus structures fllj , cap 



structures at the end of the tube |12|] and helical nanotubes . The disclinations and their 



configurations have much influence in both of their geometrical structures and electronic 



structures. [Ill] Various remarkable examples have been reported on the crucial role of the 
disclination. However, it is still unclear that what are the essential physics governing their 
effects on the electronic states of nano-structures. 

In the preceding papers, we calculated the conductances of the junctions between the 
two metallic nanotubes by using a tight binding model including only the tt orbitals. HH 
Our calculations based on Landauer's formula have revealed a remarkable scaling law; the 
conductances of the junctions are determined almost only by two parameters in a certain 
Fermi energy region. || Especially just at Fermi level of the undoped system, the conduc- 
tance is almost determined only by the ratio of the circumference of the nanotubes. || The 
result indicates that only the size and shape of the junctions are necessary to determine 
their conductances. Therefore the helicity of the honeycomb lattice is not relevant to deter- 
mine the conductance, so far as the junctions of the metallic nanotubes are concerned. It 
suggests that some continuum theory ignoring the atomic details but including only the size 
and shape of the junctions is valid, which would give some transparent view point of the 



scaling law. 

Very recently, Matsumura and Ando proposed the effective mass approximations as one 
of such continuum theories. |15| They applied it to the electron transport of the nanotube 
junctions and made the physical meaning of the scaling law somewhat clear. But their 
discussion is limited just to the case of the fixed Fermi energy case of the undoped system 
and the conductance is not given as an analytical form, as far as we know. In this paper, 
their discussion is generalized to the systems with general Fermi levels. Furthermore, the 
final result of the conductance is represented by a closed analytical form which gives a 
transparent view points. 

First of all, we explain the Bloch state of a monolayer graphite forming the single wall 
nanotube by a simple tight binding model and relate it to the effective mass theory. (16 



Fig.|l] shows the development map of the nanotube. The vector R represents the circumfer- 
ence of the tube. Two parallel lines perpendicular to R and parallel to the tube axis are 
connected with each other to form the tube. Here we use two pairs of the basic translation 
vectors {e*i,e 2 } and {e x /,e y i}, where e x r = (e\ + e 2 )/v^3 and ey = e 2 — e*i, to represent 
components of vectors on the development map. For example, the components of R in Fig.[IJ 
are represented as (R±, R 2 ) = (2,5) and (R x /,R y i) = (7-^/3/2,1.5). In this paper, we con- 
centrate our discussion to the metallic nanotube, so that only the tube of which Ri — i? 2 
is an integer multiple of three is considered. @ The four basic translation vectors have the 
same length which is about 0.25 nm and denoted by a hereafter. The wavefunction can be 
represented by the Bloch state as i/ja(q) — ex p(z(&i?i + ^2<?2))' ? /'a(0) for sublattice A and 
i/jb(q + t) — exp(i(A; 1 g 1 + fc 2 g 2 ))'?/' B (f') for sublattice B, where f = (e*i + e 2 )/3, q is the 
position of the atom with integer components q\ and g 2 . When the metallic nanotube is 
not doped, i.e., the 7r band is half filled, the Fermi energy locates at the K and K' corner 
points in the 2D Brillouin zone: corresponding wavenumbers (ki, k 2 ) are (27r/3, — 27r/3) and 
(— 27r/3, 27r/3), respectively. The corresponding energy position, i.e., the Fermi level of the 
undoped system, is taken to be zero hereafter. When the wavenumber k is near the corner 
point K, the wavenumber k' measured from the K point, (k[, k' 2 ) = (— 27r/3 + hi, 2tt/3 + fc 2 ), 



is small so that the phase factors can be linearized as exp(ikia) = w exp(ik[a) ~ w(l + ik[a) 
and exp(ik 2 a) = w' 1 exp(ik' 2 a) ~ u> _1 (l + ik' 2 a) , where w = exp(i2n/3). Then Schrodinger 
equation of a simple tight binding model for the Bloch state becomes 

Eijj A (q} = -il{wk' 2 + w~ 1 k' 1 )atljB(q + r) (1) 

and 

EipB{q + t) = i^{w" l k' 2 + wk'^aipAiq) ■ (2) 

Here 7 is the hopping integral between the nearest neighboring sites, which is about —3 eV. 
In this tight binding model, only the ir orbital is considered and mixing between the a and 
the it orbital caused by the finite curvature is neglected. By using ip' B (q) = (—i)ipB(<i + t) 
instead of ipB-, one can obtain 

^3 

EipA(q) = -7Tia(K> ~ ik'y,)il>' B (q) , (3) 



and 



y^3 
2 



W b (q) = ^a(k' x , + ik> y/ )Mq) , (4) 



where k' x , = (k^ + k^)/ v^3 and k' y , = k' 2 — k[. The solution of these equations shows the linear 
dispersion relation, 

\k'\ = ±-^- • (5) 

For the one dimensional band which intersects the K point, the periodic boundary condition 
around the circumference is R\k[ + R 2 k 2 = 0- Fro m this condition, one can show that phase 
difference between A sublattice and B sublattice is represented by 

^b(q)/Mq) = ±iexp(i0) , (6) 

where 9 is the angle of R with respect to e x > measured anti-clockwise as shown Fig.[l| . 

In the effective mass theory, the wavefunction is represented by F^ B (q)w^ qi ^ q2 \ where 
F^b is the envelop wavefunction and w <yqi ~ q2 ' > is the Bloch state wavefunction at the K point. 



When the Fermi level Ep is close to zero, the differential equation for Fj£ B is obtained from 
eq.(|3]) and eq.(|J) along a parallel discussion by replacement k' x , — > —id x i and k' y , — > —id y i as 
below. 

(-d yl - id^F* = kFf (7) 

(d y ,-id xl )FK = kFg (8) 

Here k = \k'\ and the isotropic linear dispersion relation ([|) is used. When F B = 
exp(—i9)F B < is used instead of F B , the above equations are transformed into 

(-dy-idjFjf = kFf , (9) 

and 

{d y -id x )Ff = hF B < , (10) 

where (x, y) is defined by 

d y - id x = {d y , - id x <) exp(-i0) . (11) 

The above transformation corresponds to the rotation of coordinate axes from {e x t,e y >} to 
{g x , e y } by the angle 9 anti-clockwise so that the x axis is parallel to the circumference 
vector R as shown in Figfj]. In the following discussion, the (x, y) coordinates are used 
and Fg is written simply as F B for simplicity. It can be seen from eq.(|j) that the Bloch 
state wavefunction for the one dimensional band intersecting the K point is represented by 
{F^,Fg) = (l,±i). The equations of the envelop wavefunctions F^' B for the K' corner 
point can be easily obtained in a similar way as 

(dy-id^Fg' = kFf , (12) 

and 

(-dy - id x )Ff = hFf . (13) 
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Hereafter the envelop wavefunctions F's are simply called the wavefunctions. 

The wavefunction in the nanotube is given by the plane wave; exp(i(k x x ± k y y)). When 
the x direction is taken to be parallel to the circumference of the tube, k x is quantized as 



k x (n) = 2nn/R and k y is given by k y {n) = yk 2 — k x (n) 2 . Here n is an integer representing 
a number of nodes around the circumference and R is the circumference of the tube. When 
k y (n) is a real number, the channel n is open and the corresponding wavefunction is extended, 
otherwise the channel is closed and the wavefunction shows exponential grow or decay. The 
number of the open channel is called the channel number. When the Fermi energy is zero, 
only the channel n = is open, and therefore the channel number is kept two irrespective 
of R. 

The electronic states near the Fermi energy for the undoped system (Ep = 0) govern the 
electron transport, so discussion in this paper is concentrated to this Fermi energy region 
where the channel number is kept two. In order to discuss the wavefunction in the junction 
part, the polar coordinates (r, 6) is useful. Its relation to the coordinate (x, y) is usual 
one,i.e.,r = \/x 2 + y' 2 , tan# = y/x. Fig.^ is the development map of the nanotube junction 
where the coordinate (x, y) is defined. [[7|] A heptagonal defect and a pentagonal defect are 
introduced at A(=B) and C(=D), respectively. Therefore the indexes '7' and '5' are used 
to represent the thinner tube and thicker tube, respectively. The equilateral triangles 'A 
OAB' and 'A OCD' with bases 'AB' and 'CD' have common apex 'O', which is taken to be 
the origin of the coordinate (x, y) in this paper. Then the wavefunction satisfies the wave 
equation (z 2 d 2 — zd z — <9J — z 2 )F = 0, where z = kr. The solution is represented by Bessel 
functions J m and Neumann functions N m as 

oo 

F= e tme (c m J\ m \(z) + d m N\ m \(z)) . (14) 

m=— oo 

Matsumura and Ando have found that the wavefunction should satisfy the following 



boundary conditions in the junction part. [15 



Ff(z,e + 7r/3) = -Fg(z,9) (15) 
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F§(z,0 + */3) = ~Ff (z,9) (16) 



FZ(z,e + n/3) = FX\z,6) (17) 



Fg(z,e + Tr/3) = wF%{z,e) (18) 



In these equations, w = e l27r ^ 3 . From eg . (|l5|) and eg . ([TBI) , terms in eg. (14) for F^ and F 1 



7< K 

B 



are not zero only when m = 31 + 2 (/ =integer ). Because the open channel n = has 
no node along the circumference, it is fitted to the components with smaller |m| in eg. (14) 



better than to those with larger \m\. So we assume that one can neglect all the terms except 
those with I = and I — — 1 in eg.(|l4]). Then the wavefunctions can be written as 



Ff =e 2W f 2 (z)+e^f 1 (z), (19) 



and 



f k = _jmV h{z) _ e-^h{z) , (20) 



where 



fm{z) = c m J m (z) + d m N m (z) (m=l,2). (21) 
From eg.(^|) and eg. ([13]), the other two wavefunction F§ and Ff- can be derived from 



Ff and Ff as 



Ff = -ie tf / a (z)+ie- otf /i(z) , (22) 
F* = ie^f 2 (z) - ie-^+f , (23) 



where 



/i(z) = ci J 2 (z) + diiV 2 (z), 
/ 2 (z) =c 2 J x (z) + d 2 N l (z) , 

7 



(24) 



by using well known recursion formula of the Bessel functions and Neumann functions. It 
is easily confirmed that eq.(|22|) and eq.(|23|) satisfy the boundary conditions eq.(|D]) and 
eq. (|i8|) . The amplitude of the open channel in the tube, which is denoted by a, is obtained 
from the wavefunctions in the junction part as 



A" 
7± 



V2R 



(25) 



and 



a 



K 



(26) 



5± ~ ^2R b Jc 

for the K point. The indexes + and — mean directions along which the electronic waves 
propagate as shown in Fig.^|. R$ and R 7 are the circumferences of the thicker tube and 
thinner tube, and (f) is angle between 'AB' and 'CD' in the development map FigfJ. 
Equations for the K' points are obtained from eq. ( ^5l) and eq. (|26|) by replacing i and <fi with 
—i and — 0, respectively. To simplify the calculation, the integrations in the above equations 
are transformed as fadx' — *> R$ j_£ _,d9 and dx — > R7 S_I d9. If the wavefunction 

3 7T (p S TT 

in the junction varies slowly near r = R§ and r = R 7 along the radial directions, this 
replacement can be allowed. 

The relation among the amplitudes of the open channel in the thicker tube, c?5 = 



5L), those in the thinner tube, a 7 = 1 (otf + , af+, af_, af^_) and the wave- 



functions in the junction part, c = * (c 2 , d 2 , Ci, d\) is summarized in the followings. 



0C5 



/ i? 5 n(0)MA(0)L(fci? 5 )c 



(27) 



a 7 - 

where M are constant matrixes as 

/ 

M = 



'R 7 P 7 ML(kR 7 )c 



w ~^§iw^ 

2 

% 

-10 -^-i 



(28) 



(29) 



S 



A(0) is a diagonal matrix, where Ax,i = Ag 3 = exp(— and A 2)2 = A^ 4 = exp(— 2i(p). Pj 



and Ps((j)) are defined by eq.fl2"5|) and eq. (p6|) , respectively. The matrix elements of L(z) are 
Ln = L 33 = Jx(z),L 12 = L 34 = Ni(z),L 2 i = L43 = J 2 (z) and L 22 = L M = N 2 (z). The other 
matrix elements of L(z) are zero. The vector a-j is given by a linear transformation of c?5 
defined by eq.(P?D and eq.(^Sj). The transformation matrix T (i.e., a-j = T0J5) includes three 
independent parameters, i.e., (3 = R7/R5, z$ = kR 5 and 0. The parameter z 5 is related to 
the Fermi Energy Ep as follows. When \k\ is near zero, channel number is always two which 
is independent of the radius of the nanotubes. But as \k\ increases, the channel number 
increases firstly in the thicker tube, when \k\ exceeds k c = 2tt/R^ . Owing to the linear 
dispersion relation eq.(|5]), z 5 = 2nk/k c = 2nEp/E c holds where E c is the threshold Fermi 
energy corresponding to k c . The transmission rates are calculated from the transfer matrix 
T, and the conductance a is obtained by Landauer's formula as 

24 



a 



{/3£Li£? =1 (3/4)^- 2 X^} + 6 



(30) 



where 



Xij = nztiJiiPz^Njizn) - NiiPz^Jjizs)} . (31) 

The obtained conductance has a remarkable feature that it does not depend on the angle 
<p. It is consistent with the scaling law with the two parameters Ep/E c and (3 = -R7/-R5 in 
Ref. 0. When Ep = 0, a has a simpler form as 

a = 8/((3 3 + (3~ 3 + 2) . (32) 

Eq. (|52"D is also fitted well to the numerical data in Ref. which shows one parameter scaling 
law. 

Fig.^ and Fig.[| show the comparison between the conductances obtained by the tight 
binding model and those by eq.(|30"D for the (2i, 2i + 3)-(10,13) junctions, where the same 
notation as that of Ref. is used to specify the junctions. The horizontal axis is l-E^I 
with the unit of the absolute value of the hopping integral |7|. The figures are shown in the 



energy range where the channel number is kept two, i.e., \Ep\ < E c . Agreement between the 
two is quite good. But difference is larger for positive Fermi energy than that for negative 
Fermi energies. This difference is explained as the effect of discreteness of the lattice as 
follows. The wavefunction tends to have opposite signs between neighboring sites when the 
Fermi energy is positive. But such alternating sign arrangements are not possible for the 
paths with odd number steps surrounding the odd membered ring defects. Such effects come 
from the discreteness of the lattice and can not be treated by the continuum model assumed 
in the the effective mass theory. It should be also noted that a drop of a near Ep = E c 
appears in the tight binding model, while it is not reproduced by eq . (|30|) . It comes from 
the effects of the closed channel which are not included in eq. (|30|) . As Ep approaches E c , 
the decay length of the closed channel n = 1 increases and becomes more extended so that 
it has more important role in the transport. When the values of (3 = R7/R5 approaches to 
unity 1, the conductance does not depend on Ep taking the values close to the maximum 
value. As the parameter f3 decreases, the conductance a near Ep = decreases and peak 
structure appears at Ep somewhat lower than E c . When f3 increases further, the height of 
the peak is lowered. The value of z 5 corresponding to the Fermi energy at the peak is about 
half of the first period of J^{z), and 3/4 of that of J\{z) so that we speculate that the peak 
structures are formed by resonance of these radial wavef unctions. More detailed discussions 
will appear elsewhere. 

In this paper, the scaling law of the conductance of the nanotube junctions previously 
found by the tight binding model can be reproduced and analyzed by the effective mass 
theory. It provides the analytical representation of the conductances and a clear physical 
interpretations. It can be also applied to various structures relating to the nanotubes. 

We would like to thank H. Matsumura and T. Ando for their useful suggestions. This 
work has been supported by the Core Research for Evolutional Science and Technology 
(CREST) of the Japan Science and Technology Corporation(JST). 
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FIGURES 



FIG. 1. Development map of the nanotube. 

FIG. 2. Development map of the nanotube junctions. It is similar to that of Ref. 0. Lines 'EB' 
and 'FA' are parallel. Lines 'DG' and 'CH' are also parallel and their angle with respect to 'EB' 
and 'FA' is denoted by cp. The lines 'EB','BD' and 'DG' are connected with the lines 'FA', 'AC 
and 'CH', respectively. The rectangles 'EBAF' and 'DGHC form the thinner tube and the thicker 
tube, respectively. 'BD' is the rotated 'AC by angle of 60 degree and the quadrilateral 'ABDC 
forms a junction part with a shape of a part of a cone. A heptagonal defect and a pentagonal 
defect are introduced at A(=B) and C(=D), respectively. 

FIG. 3. The conductances of the junction between (2i, 2i+3) tube and (10,13) tube for negative 
Fermi energy. Corresponding values of i are attached to each graphs. The horizontal axes are the 
absolute value of the Fermi energy in units of the absolute value of the hopping integral \ j\ = —7. 
The hopping integral 7 is taken to be negative. They are calculated by eq. (|30|) , and the tight 
binding model, each of which are represented by dotted lines and plots connected with solid lines. 

FIG. 4. The same figure as Fig. 2 for the positive Fermi energy. 
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